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We explore the coherent destruction of tunneling (CDT) in a lattice array under selective in- 
phase harmonic modulations, in which some selected lattice sites are driven by in-phase harmonic 
oscillating fields and other lattice sites are undriven. Due to the occurrence of CDT, if the driving 
amplitude A and the driving frequency uj are tuned to satisfy the zeroth-order Bessel function 
Jo(A/lu) = 0, the driven lattice sites are approximately decoupled with the undriven lattice sites. 
The CDT even takes place in lattice systems with high-order couplings between non-nearest lattice 
sites. By using the CDT induced by selective in-phase harmonic modulations, we propose a scheme 
for realizing directed transport of a single particle. It is possible to observe the CDT in the engineered 
optical waveguide array, which provides a new opportunity for controlling light propagation and 
designing switch-like couplers. 

PACS numbers: 42.65.Wi, 42.25.Hz 



I. INTRODUCTION 

Controlling quantum states via external driving fields 
is a very popular approach for quantum manipulation, 
which is a basic tool for implementing quantum technol- 
ogy. The external driving fields have been used to ma- 
nipulate various systems including spins, quantum dots, 
cold atoms and Josephson junctions etc [l], The co- 
herent destruction of tunneling (CDT) Q and dynamical 
localization (DL) [4j are two seminal results in this field. 
Moreover, some important generalizations of CDT, such 
as selective coherent destruction of tunneling (SCDT) Q, 
non-degenerate CDT Q, nonlinear CDT Q, and many- 
body CDT [8J, have also been explored. Single-particle 
CDT and DL which have been observed exper- 

imentally in optical traps, attract considerable atten- 
tion due to their potential applications in quantum mo- 
tor [ill . Il2j and quantum information processing fl3l j . It 
has also demonstrated a directed transport in a driven 
bipartite lattice based on the idea of CDT 14 1. 

Engineered optical devices offer a new platform for 
studying quantum dynamical control [lj|. In contrast 
with a quantum system, an optical device provides the 
unique advantage of a direct mapping of the temporal 
evolution of the quantum system into spatial propagation 
of light waves in engineered lattices. DL in waveguides 
and CDT in directional couplers are visualized either by 
periodically bending waveguides [H, [TtJ , or by an out- 
of-phase harmonic modulation of refractive index [l8j . 
These works focus on considering photonic structures 
with out-of-phase modulations of refractive index of ad- 
jacent waveguides [l8| . For a single-band quantum lat- 
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tice system of periodic modulations on left part of the 
lattice [l9j . it has demonstrated the occurrence of dissi- 
pationless directed transport. However, the question of 
whether quantum tunneling can be suppressed or even 
inhibited via selective in-phase modulations is still not 
clear. Here, the selective in-phase modulation means 
that the refractive indices for the waveguides oscillate 
in-phase in selected waveguides, while they are propa- 
gation invariant for all other waveguides. Within the 
context of a driven optical lattice, the in-phase modula- 
tion is termed when the selected lattice sites are shaken 
up and down in-phase while all other lattice sites re- 
main static. In addition, a major limitation of previ- 
ous schemes for quantum dynamical control is that it is 
only suitable for systems under nearest-neighbor tight- 
bing approximation (NNTB), but fails to explain lattices 
with long-range interaction. 

In this article, we analyze how CDT is induced by se- 
lective in-phase harmonic modulations and discuss some 
possible applications of CDT in quantum motor and 
switch-like coupler. Due to CDT occurs only between the 
modulated lattice site and the unmodulated ones, one can 
get the precise location of a particle in the lattice arrays. 
Under the selective in-phase harmonic modulations, the 
coupling between a driven lattice site and an undriven 
lattice site is rescaled by a factor of the zeroth-order 
Bessel function, while the coupling among the driven lat- 
tice sites or the undriven lattice sites keeps unchanged. 
By tuning the driving frequency and the driving ampli- 
tude, the effective tunneling between a driven lattice site 
and an undriven lattice site can be completely suppressed 
and therefore the lattice array divides into some sets of 
disconnected lattice. By using CDT, it is possible to se- 
lectively move the quantum particles and design a new 
type of optical switch-like devices. 

The article structure is as follows. In the next section, 
we give the Hamiltonian for the quantum lattice systems 
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studied in this work. In section III, we explore how CDT 
is induced by the selective in-phase harmonic modula- 
tions. We obtain some analytical results which are con- 
firmed by our numerical simulation. In section IV, we 
discuss some possible applications of CDT in quantum 
motor and switch-like coupler. In the last section, we 
briefly conclude our results. 



II. HAMILTONIAN 

We consider a single particle in an array of lattice sites, 
which is described by a single-band tight-binding model. 
In addition to the tunneling between nearest-neighboring 
sites, we take account into the high-order coupling be- 
tween non-nearest-neighboring sites. The Hamiltonian 
reads as 

H = Yl + h - c -) + £ v ^)+2 h i + h - c -) 

3 j 

+ X>i(*)a]aj, (1) 
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where S7 and v are the tunneling amplitudes connecting 
nearest-neighboring sites and next-nearest-neighboring 
sites, respectively. The (dj) are the particle creation 
(annihilation) operators in the j-th site. 

The in-phase harmonic oscillating fields are applied to 
modulate the on-site energy £j(t). We consider a type of 
selective in-phase harmonic modulations of £j(t) = for 
the undriven lattice sites and £j(t) — As'm(ujt) for the 
driven lattice sites. Here, A is the driving strength and 
lu is the driving frequency. The case of £i{t) — £i+i(t) = 
Asm(ujt) and £j(t) = for j ^ i, z+1 is a typical example 
of selective in-phase harmonic modulation, where the i- 
th and (i + l)-th sites are modulated in phase while all 
other sites are unmodulated. 



III. COHERENT DESTRUCTION OF 
TUNNELING INDUCED BY SELECTIVE 
IN-PHASE MODULATIONS 

Taking the Wannier state \j) — at |0) localized in the 
jth site as basis, we expand the quantum state ^(t)) of 
system (fTJ in the form 



!*(*)> = £»i(*)|j>. 



(2) 



where dj(t) represents the occupation probability am- 
plitudes in the site j, with the normalization condition 

E,-KI 2 = i- 

From the Scrodinger equation, idt\^f(t)) — H\^(t)), 
the evolution equation for the probability amplitudes 
a,j(t) reads 

= ^ a i- 1 + a i+ 1 ) + u ( a j'-2 + a j+2) + Sj(t)aj. (3) 



The transformation a,j = bj exp[— i J Sj(t)dt] yields 

db f 
i-^j- — fi6j_icxp[— i I (sj-i — £j)dt] 

+£lbj+i exp[— i J {sj+i — £j)dt] 



+vbj + 2 exp[— i / (ej+2 — £j)dt] 



-vbj-2 exp[ 



(£j_ 2 -£j)dt}. (4) 



For a iV-site system under in-phase modulation, the 

values of £j{t) (j = 1, , N) are only either £j(t) = 

or £j(t) = Asm(u>t), so the values of £j{t) — £j>(t) for 
j j' should be either or ±^4sin(o;i). The former 
case of £j (t) — Eji (t) — indicates that the j-th andj'-th 
sites are simultaneously modulated in-phase or unmodu- 
lated, and the latter case of Sj (t) ~ £y (t) = ±A sin(wt) 
means that one site is modulated while the other is un- 
modulated. For the case of £j(t) — £j'(t) = 0, one 
finds that the effective tunneling parameters between 
the jth and j'th sites remain unchanged. For the case 



Of £j (t) 



'(*) 



±Asin(wt), using the expansion 
exp[±iAcos(ujt)/uj] — ^ fe (±i) fe Jfe(A/w) exp(±ifcwt) in 
terms of Bessel functions and neglecting all orders ex- 
cept k = in the high frequency region, one finds that 
the effective tunneling parameters between the j-th and 
j'-th sites are renormalized by a factor of Jq(A/u>). 

In our analysis, we focus on the dynamics for three 
sites. However, similar behaviors may appear in other 
lattice systems as well. For a three-site system, the cou- 
pled equation §5§ reads as 

da\ 

i—— = ila 2 + va 3 + £i(t)ai, 
dt 



.da,2 



i—— = flai + Qa 3 + £2(^)02, 
dt 

da,3 . . 

i—j- = Ma,2 + va\ + £ 3 (tja 3 . 



(5) 



We consider two typical in-phase modulations for the 
above three-site system: (a) £\{t) — Asin(uJt), £2{t) — 
£ 3 (t) = 0; and (b) £ X {t) = 0, £ 2 {t) = £ 3 (t) = Asm(wt). 
In the high-frequency limit, taking time average of high 
frequency terms, we find the same effective equation of 
motion for both modulations, 



dbx 
'~dt 

1 dt 

' dt 



nj (A/uj)b2 + vJ Q (A/uj)b 3 , 
Jo (A/u) h + Slba, 
£lb 2 + vJ (A/u)bi. 



(6) 



According to the Floquet theorem, the approxi- 
mate Floquet solutions of Eq.© can be constructed 
as dj = exp[—if£j(t)dt]bj(t) (j = 1,2,3) with 
{h(t),b 2 (t), b 3 (t)) T - {B 1: B 2 , B 3 ) T exp{-iEt) being the 
stationary solution for Eq. ([B]). 
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If Jq(A/uj) = 0, bi(t) does not couple with b 2 (t) and 
6 3 (t). This means that the population in state |1) does 
not change with time. However, and bs(t) are still 
coupled and so that the populations in states |2) and |3) 
may exchange. There are several roots for Jq{A/lu) = 
such as A/u ~ 2.4 and 5.5. For a system of type-(a) 
modulation and Jq(A/lo) = 0, the approximate Floquet 
solutions for Eq. (|S|) read as 



\ui(t)) = exp[iAcos(ujt)/uj]\l] 
1 ,_, 1 



M*)) 
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[0|1> 
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V2 
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+ 0|2)+0|3>, 
|3)] exp(— iilt), 

|3)] exp(ifii). (7) 



Clearly, the Floquet state \ui(t)) is localized in the site 1, 
and 1^2,3(4)) are localized in the sites 2 and 3 with equal 
probability. In Fig. 1, we show the dependence of the 
quasienergy on the ratio A/uj. The quasienergy spectra 
are identical for type- (a) (solid lines) and type-(b) (cir- 
cles) modulations. The vertical dashed line corresponds 
to J (A/u) = 0. In Fig. [1] (b), for the Flqouet state 
at Jo(A/ui) = 0, we show the probability of finding the 
particle in different lattice sites. Generally, the Floquet 
states are either localized in the site 1 or equally dis- 
tributed on both site 2 and site 3, which confirm the 
approximate Floquet solutions very well. 
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FIG. 1: (color online) (a) Quasienergies versus A/uj. The 
vertical dashed line labels the position of Jo(A/u>) = 0. The 
solid lines and circles correspond to type-(a) and type-(b) 
modulations, respectively, (b) The probability distribution 
of the Floquet states for the system of Jo(A/uj) — labeled 
by the vertical dashed line in (a). The other parameters are 
chosen as £1 = 1, v = 0.2 and to — 10. 

Our numerical simulation is implemented by directly 
integrating the time-dependent Scrodinger equation ([5]) 
with type-(a) or type-(b) modulations. For different 
initial conditions, we calculate the probabilities Pj = 
l a il 2 '0 — 1)2,3) of finding the particle in lattice site 
j. In Fig. [Ha), we show the time evolution for the 
system of A/lu ~ 2.405 and the particle is initially lo- 
calized in the site 1. Since the tunnelling between the 



site 1 and the other two sites are suppressed, the par- 
ticle is frozen in the site 1. On the other hand, if we 
start a particle in the site 2, the particle indeed simply 
tunnels back and forth between the sites 2 and 3, see 
Fig. [2](b). These results indicate that only the tunnel- 
ing between the site 2 and 3 is allowed. This is a kind of 
selective CDT, which is also evidenced by the probabil- 
ity distribution of the Flqouet states shown in Fig. [1] (b) . 
One can use the effective Scrodinger equation ([6]) to ex- 
plain the SCDT. Under the condition of Jq(A/u) = 
and b\ = 1, b% = 63 = 0, one can immediately give the 
analytical solution as 61 = l,b% = ba = 0. Similarly, 
if the initial condition takes b\ = 0,62 = 1;&3 = 0, the 
analytical solution is given as b\ — 0, 62 = cos(Oi) and 
&3 = —i sin(ilt). 




FIG. 2: (color online) Time evolution of the probability dis- 
tribution for the system of A/uj ~ 2.405. The particle is 
initialized in site 1 [column (a)], and site 2 [column (b)]. The 
modulations and the parameters are the same as the ones for 
Fig.ffl 

Although the above analysis is applied to three-site 
systems, similar analysis can be applied to systems of 
arbitrary number of lattice sites and similar CDT phe- 
nomena can be found. The effective inter-site tunnelings 
are illustrated by the cartoon shown in Fig. [3] If system is 
driven by a high-frequency field such that Jq(A/qj) = 0, 
as illustrated in Fig. [31 then the coupling parameters 
between the unmodulated and modulated site are ef- 
fectively rescaled by the factor J , while the couplings 
among modulated sites or unmodulated sites keep un- 
changed. The upper cartoon in Fig. [3] represents the 
systems of modulation on site i only, in which the site i 
becomes a true trap for the particle, since tunneling in 
both directions are suppressed, which is marked by the 
red cross. If in-phase modulations are applied to both 
sites i and i + 1, the lattice array divides into three sets 
of disconnect sites (...,« — 3, i — 2,i— 1), + 1) and 
(i + 2, i + 3, ...), see the lower cartoon. If the particle is 
initialized in site i or i + 1, Rabi oscillation will occur 
between site i and i + 1, and the frequency of this oscilla- 



tion is determined by the value of the original tunneling 
strength CI. 
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FIG. 3: (color online) The schematic diagram for the dy- 
namics of lattice arrays under selective in-phase harmonic 
modulations. Upper cartoon: the modulation is only ap- 
plied to site i (marked by £), that is, £i(t) — Asin(u)t) and 
~ for (j 7^ i). Lower cartoon: in-phase modula- 
tions are applied to sites i and i + 1 (marked by £), that 
is, Ei{t) — Ei+i{t) = Asm(ujt) and £j(t) = for (j ^ + 1). 



IV. POSSIBLE APPLICATIONS 

Based on the CDT in a lattice system under selective 
in-phase modulations, we propose a scheme for directed 
single-particle transport, which can be used to realize 
a quantum motor. The transport process is illustrated 
in Fig. |U in which A/uj ~ 2.405 and w/fi = 10. The 
particle is initially loaded into site i and in-phase mod- 
ulations are applied to sites i and i + 1, see Fig. [H The 
Rabi oscillation between sites i and i + 1 has a period 
t = Tr /CI. At t = t/2, the particle arrives at site i + 1 
and the in-phase modulations are suddenly changed to 
address sites i + 1 and i + 2. Therefore the particle enters 
into a new Rabi oscillation between sites i + 1 and i + 2 
with the same period r. Similarly, at t = r, the modu- 
lations are suddenly switched to address sites i + 2 and 
i + 3, the particle has completely tunneled to site i + 3 
after another time interval of r/2. If the in-phase modu- 
lations are repeatedly applied to the lattice site occupied 
by the particle and its right-side neighbor at t = nr/2 for 
n = 0,1,2,3..., the particle will propagate solely to the 
right. Similarly, the in-phase modulations are repeatedly 
applied to the lattice site occupied by the particle and 
its left-side neighbor at t — nr/2 for n = 0, 1, 2, 3..., the 
particle will propagate solely to the left. 

The directed quantum transport based upon CDT pro- 
vides an optional scheme for designing a quantum motor. 
To illustrate the possibility of realizing a quantum motor 
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FIG. 4: (color online) Schematic diagram for transporting 
a particle from left to right in a lattice array. The symbol 
£ denotes the high-frequency modulation e = Asin(ujt), in 
which A/uu ~ 2.4. After a time interval of t = r/2 with 
r = n/Q, the particle has completely transported from one 
site to its right-side neighbor. 

via the directed quantum transport, we show our numeri- 
cal results of a single particle in a 11-site system with the 
procedure described above, see Fig. [5l Obviously, the di- 
rected single-particle transport is perfectly confirmed by 
our results. 
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FIG. 5: (color online) Simulation of a single particle in a 11- 
site system. The situation is illustrated in FigU 

In addition to the application in designing a quan- 
tum motor, the CDT may be used to realize a quan- 
tum beam splitter. If the in-phase modulations are 
initially applied to the sites i — 1, i and i + 1, that 
is, £,;_i = £j = Ei + \ = Asu\(iot) and Ej(t) = for 
j ^ i — 1, i, i + 1. Therefore the three sites i — 1, i and 
i + 1 are disconnected with other sites, and the couplings 
among the three sites keeps unchanged. The dynamics 
can be well effectively described by Eq. ([S]) for the case 
of no modulations (i.e., E\ = £2 = £3 = 0). Initializ- 
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ing the system into a state of 02 = 1 and a\ — 03 = 0, 
we can derive the analytical solution for Eq. with 
e% = £2 = £3 =0, which has equal probability ampli- 
tudes for the first and third sites(Pi = P3). The particle 
has same probability for tunneling to sites i—1 and i + 1 
and therefore the state splits into two parts. At the time 
of zero probability for finding the particle in the site i, 
we switch off the in-phase modulations on the site i and 
simultaneously switch on in-phase modulations on sites 
i — 2 and i + 2. After a time interval of r/2, one will find 
the particle in site i— 2 or i + 2 with equal probability. Af- 
ter a time interval of t/2, we switch off the modulations 
on sites i — 2 and i + 2 and simultaneously switch on the 
modulations on sites i— 3 and i+3. Again and again, after 
each time interval of t/2, we switch off the modulations 
on sites i — k and i + k and simultaneously switch off the 
modulations on sites i— (fc+1) and i + (k + l). Under this 
sequence of applying the in-phase modulations to some 
certain lattice sites, the particle initially localized in the 
central site will form an equal-probability superposition 
of left and right propagating components. This means 
that a perfect quantum beam splitter is realized by the 
CDT caused by the in-phase modulations. 
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FIG. 6: (color online) A quantum beam splitter for a particle 
in a 11-site system. If the in-phase modulations initially acts 
on the site 6 and its nearest neighbors (the sites 5 and 7) , the 
state splits into two equal parts propagating along opposite 
directions. 

Although some other good schemes for directed trans- 
port of single-particle has been previously proposed in a 
driven lattice [14|, we stress here that one advantage of 
our approach is that it can be applied to lattices with 
high-order coupling among non- nearest sites. 



V. OPTICAL REALIZATION 

By mapping the temporal evolution of quantum sys- 
tems into the spatial propagations of light waves, the 



engineered waveguides offer an alternative platform to 
investigate the classical wave analogues of a wide variety 
of coherent quantum effects. Below, we discuss how to 
simulate the CDT induced by selective in-phase modula- 
tions via the engineered waveguides. 

The propagation of a cw light wave along the z axis is 
described as the Schrodinger equation for its dimension- 
less field amplitude 6(x. z)[18L l20j 



(8) 



Here (x and z) are the normalized (transverse and longi- 
tudinal) coordinates, and p describes the refractive index 
contrast of the individual waveguide. The refractive in- 
dex distribution of the waveguide lattice is given by 



A I 



R(x,z) = [l + fj(z)]exp 

j = -M 

fj( z ) = Fjfism^z), 



x - jw s 



w. 



■ (9) 



with the waveguide spacing w Sl the channel width 
w x , the longitudinal modulation amplitude n, and the 
modulation frequency u). The super-Gaussian function 
exp(— x 6 /w®) describes the profile of individual waveg- 
uides with widths w x . 2M + 1 is the total number of the 
waveguides. The function Fj determines the modulation 
type. Here we consider the selective in-phase harmonic 
modulations, such that the value of Fj can only be either 
Fj = 1 or Fj = 0. The case of Fj = 1 indicates that the 
j-th waveguide is modulated, while Fj = means that 
the j-th waveguide is unmodulated. 

Under the tight-binding condition, one can expand the 
total electric envelope in form of 



(x,z) = aj(f>j(x) exp(-i0 o z) 

A) 

#0 



(10) 



(j)j(x)H (j)j(x)dx, 

2 » 
2c5a? 



1 d' x ^ 

j=-M 



x - ]W S 



where 4>j{ x ) are the eigen-modes for the j-th waveguide, 
and the expanding coefficients are normalized to one, 
12j \ a j\ 2 = 1- Substituting the expansion (10) into the 
Schrodinger equation ((8|) and taking into account cou- 
plings between next-nearest and nearest waveguides, one 
can obtain the resulting coupled-mode equation in the 
same form as Eq. with t replaced by z. Therefore, 
the temporal evolution of the probability amplitudes dj 
obeying Eq. © can be simulated by the light propaga- 
tion in an array of tight-binding waveguides with longi- 
tudinal modulation, in which the temporal evolution of 
a,j is mapped into the spatial evolution of the modal field 
amplitudes of light waves in the j-th waveguide along the 
axial direction x. The light power distribution P(x, z) 
corresponds to the probability distribution |a,, (i)| 2 . 
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To demonstrate the quantum control based on the 
CDT, we simulate the three-channel couplers [M = 1) 
by integrating the field propagation equation dHJ. In our 
numerical simulation, the initial states are chosen as the 
lowest Wannier modes for isolated individual waveguides, 
and the parameters are set as w x = 0.3, w s = 3.2 and 
p = 2.78. As same as the units used in Refs. [H, HH, w x 
and w s are in units of 10 /im, and p = 2.78 corresponds 
to a real refractive index of 3.1 x 10~ 4 . 

According to the analytical analysis in Section III, if 
not all three waveguides are modulated simultaneously 
and modulations satisfy J (A/uj) = 0, the three-channel 
couplers can be decoupled into a coupled two- waveguide 
system and an isolated single- waveguide system. To com- 
pare with the sub-systems decoupled from the original 
system, we also simulate the tunneling dynamics of light 
in an unmodulated two-waveguide system corresponding 
to the two waveguides decoupled from the original sys- 
tem, see Figs. [7[b). 

Similar to the discrete systems, the in-phase modu- 
lations applied to the continuous system @ may sup- 
press the tunneling between modulated and unmodulated 
waveguides for some specific modulation amplitudes and 
frequencies. We find that the tunneling between modu- 
lated and unmodulated waveguides is strongly suppressed 
if the modulation frequency u = 3.45 x (27r/100) and 
/.t = 0.4. This strong suppression of the tunneling be- 
tween modulated and unmodulated waveguides is a new 
extension of CDT. 

We simulate a three-channel coupler under two typical 
in-phase modulations: (i) in-phase modulation is only 
applied to one waveguide [see Figs. [7] (c) and (e)], i.e. 
Fi = I and F-i = Fq = 0; and (ii) in-phase modulations 
are applied to two next-nearest-neighboring waveguides 
[see Figs. M (a) and (c)], i.e. F x = F_ x = 1 and F = 
0. We numerically explore the light propagation in a 
three-channel coupler under these two modulations. In 
these two figures, the left columns show the refractive 
index distributions R(x, z) and the right columns show 
the light propagation \(f)(x, z)\ 2 which are obtained by 
numerically integrating the Schrodinger equation ([8]). In 
our simulation, the parameters are set as w x = 0.3, w s = 
3.2, p = 2.78, n = 0.4 and ui = 3.45 x (2tt/100). 

In Fig. 7 (c)-(f), for a three-channel coupler under 
type-(i) modulation, we show the refractive index dis- 
tribution R(x 7 z) [see panels (c) and (e)] and the light 
propagation \<j)(x, z)\ 2 [see panels (d) and (f)] for different 
input beams marked by The first row shows R(x, z) — 
exp[— x 6 /w®] + exp[— (x + w s ) 6 /w x ] and \4>(x,z)\ 2 for the 
unmodulated two- waveguide system corresponding to the 
subsystem decoupled from the original three-channel sys- 
tem, see Figs. 7 (a) and (b). If the input beam is cen- 
tered in the middle waveguide marked by ► [see panel 
(c)], which is one of two unmodulated waveguides, the 
light periodically tunnels between the two unmodulated 
waveguides. The tunneling dynamics is almost the same 
as the one in the unmodulated two-waveguide coupler 
with R(x,z) = exp[— x 6 /w x ] + cxp[— (x + w s ) 6 /w®], see 



Fig. [T^b) . While the input beam is centered in the sole 
modulated waveguide marked by ► in Fig.[7](e), the light 
propagates in this modulated waveguide almost with- 
out tunneling into the two unmodulated waveguides, see 
Fig. [3(f). The light propagation shown in Figs. [7] (d) and 
(f) clearly shows the strong suppression of the coupling 
between modulated and unmodulated waveguides. 
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FIG. 7: (color online). Light propagation in a three-channel 
coupler under type-(i) modulation. First row: (a) the refrac- 
tive index distribution R(x,z) and (b) the light propagation 
|0(x,z)j 2 for the unmodulated two-waveguide system corre- 
sponding to the two unmodulated waveguides decoupled from 
the original three-channel system. Second row: (c) R(x, z) 
and (d) \(f)(x, z)\ 2 with the input beam centered in the middle 
waveguide marked by ► in (c). Third row: (e) R(x, z) and (f) 
\<f>(x,z)\ 2 with the input beam centered in the top waveguide 
marked by ► in (e). 



For a three-channel coupler under type-(ii) modula- 
tion, in which the top and bottom waveguides are mod- 
ulated, we show the corresponding light propagation for 
different input beams in Figs. [5] (a) and (c). Similarly, if 
the input beam is centered in one of the two modulated 
waveguides, the light only oscillates between the two 
modulated waveguides, see Fig. |8](a) and (b). While the 
input beam is centered in the middle waveguide marked 
by ► in Fig.|S](c), the light almost perfectly keeps prop- 
agating in this waveguide, see Fig. [5] (d). Obviously, 
the numerical results indicate the decoupling between the 
modulated and unmodulated waveguides. 

For a quantum lattice system, it is generally thought 
that the nearest-neighboring tunneling is more signif- 
icant than the next-nearest-neighboring one. How- 
ever, due to the in-phase periodic modulations, even 
the nearest-neighboring tunneling strength is stronger 
than the next-nearest-neighboring one, the appearance 
of significant tunneling between next-nearest-neighboring 
waveguides and almost no tunneling between nearest- 
neighboring waveguides updates this conventional under- 
standing. This provides an optional approach for directly 
transporting a particle to its non-nearest-neighboring 
sites. 



7 




500 1000 1500 500 1000 1500 



z z 

FIG. 8: (color online) Light propagation in a three-channel 
coupler under type-(ii) modulation. First row: (a) R(x, z) and 
(b) \(j>(x, z)\ 2 with the input beam centered in the top waveg- 
uide marked by ► in (a). Second row: (c) R(x,z) and (d) 
\4>{x, z)\ 2 with the input beam centered in the middle waveg- 
uide marked by ► in (c). 

For different input beams, we have also numerically 
explored the light propagation \<f>(x, z)\ 2 in a three- 
channel coupler under in-phase modulations applied to 
two nearest-neighboring waveguides. Again, the light 
propagation (not shown here) clearly indicates that the 
modulated and unmodulated waveguides are almost de- 
coupled each other and there is almost no tunneling be- 
tween them. 

Moreover, we have simulated the multi- waveguide sys- 
tems of other numbers of waveguides by directly inte- 
grating the Schrodinger equation (|5]l. The predictions 
from the coupled-mode equation ([3]) are demonstrated 
and confirmed perfectly in this optical system (not shown 
here). Therefore, such a waveguide array under in-phase 
modulations may be used as a switch-like device, which 
can function as input-dependent single-mode fibers, dual- 
channel couplers and multi-channel couplers as well. 

VI. SUMMARY 

In conclusion, we have studied how to control the tun- 
neling dynamics in one-dimensional quantum lattice sys- 



tems via CDT induced by selective in-phase modulations, 
in which only some selected lattice sites but not all lattice 
sites are driven by periodic in-phase fields. Tuning the 
frequency and amplitude of the selective in-phase modu- 
lations to some specific values, it is possible to find that 
the CDT only occurs between modulated and unmodu- 
lated sites while the tunneling dynamics among modu- 
lated or unmodulated sites keeps unchanged. In partic- 
ular, by applying in-phase modulations to next-nearest- 
neighboring sites, it is possible to switch off the tunnel- 
ing between the nearest-neighboring sites. Therefore, one 
can utilize the in-phase modulations to produce ratchet- 
like motions and even to transport directly a particle be- 
tween non-nearest-neighboring lattice sites with the high- 
order coupling among those lattice sites. 

Furthermore, we have shown that our results on CDT 
can be experimentally tested with an array of optical 
waveguides, in which the spatial propagation of light in 
an optical material system corresponds to the temporal 
evolution of state in a quantum system. Comparing with 
previous works on light propagation in optical waveguides 
under out-phase modulations [HI, [2(| , the in-phase mod- 
ulations provide a new possibility for controlling light 
propagation in periodic waveguide arrays. 
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